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Abstract 

Linear gyro-kinetic stability calculations predict unstable toroidal Ion Temperature Gradient 



^ \ modes with normalised poloidal wave vectors well above one (k$Pi > 1) for standard parameters 



and with adiabatic electrons. These modes have a maximum amplitude at a poloidal angle 9 that 

is shifted away from the low field side (9 ^ 0). The physical mechanism is clarified through the 

K^ >j use of a fluid model. It is shown that the shift of the mode away from the low field side (9 ^ 0) 

, &j reduces the effective drift frequency, and allows for the instability to develop. Numerical tests 

7— ( ■ using the gyro-kinetic model confirm this physical mechanism. It is furthermore shown that modes 

CN . with 9^0 can be important also for kgpi < 1 close to the threshold of the ITG. In fact, modes 
ON- 

I with 9^0 can exist for normalised temperature gradient lengths below the threshold of the ITG 
obtained for 9 = 0. 
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I. INTRODUCTION 



The growth rate of the ion temperature gradient mode (ITG) as a function of the nor- 
malised poloidal wave vector kgpi, has been reported many times in the literature, see for 
instance as a bell shaped curve with a single maximum. Here, kg is the poloidal com- 



ponent of the wave vector and pi = rriiVthi/ ZeB = y/1m(Til 'ZeB is the ion Larmor radius, 
with rrii the ion mass, Vtu the ion thermal velocity, Z the charge number, e the elementary 
charge, B the magnetic field strength, and Tj the ion temperature. 

In this paper we report on collisionless Ion Temperature Gradient (ITG) instabilities 
with kg pi > 1 and adiabatic electrons. It will be shown that these instabilities can exist for 
relevant Tokamak parameters. The physical mechanism of these instabilities will be shown 
to be related to the reduction of the effective drift frequency through the shift of the mode 
away from the low field side position. This mechanism is different from previously reported 
{2-6] instabilities with kePi > 1, which are unstable only in a slab or for weak toroidicity. 

This paper is structured as follows: Section II introduces the high kg ITG through nu- 
merical simulation based on the gyro-kinetic model, Section III discusses the physics of the 
instability through the use of a simple fluid model. Section IV discusses the relation with 
previously published work and, finally Section V gives the conclusions. 

II. HIGH kg ITG 

An example of a calculation using the ballooning transform [?| (the use of this transform 
is referred to below as the 'spectral case'), obtained with the gyrokinetic code GKW [8], is 
given in Fig. 1 by the dash-dotted curve that has a single maximum. The parameters of 
this, and all other simulations in this paper, are those of the Waltz standard case ion 
temperature gradient lengths R/L Ti = 9.0 density gradient length R/L ni = 3.0, electron 
and ion temperature T e = Tj, safety factor q = 2, magnetic shear s — 1, and inverse aspect 
ratio e = 0.166. The simulations use circular geometry retaining finite e effects, and the flux 
tube approximation is always applied. 

However, GKW simulations with the radial direction described using finite differences 
(simulations that use finite difference in the radial direction are referred to below as the 
'non-spectral case') show a surprisingly different behaviour for kgpi > 0.6, displaying a 
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spectrum with two maxima and having unstable modes with kgpi well above one, as is 
shown by the full line of Fig. 1. 
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FIG. 1: (Colour on-line) Growth rate (7 in units Vthi/Ro where Rq is the major radius of the 
magnetic clXls) clS cl function of kgpi. The dash-dotted (blue) line is the spectral case with #0 = 0, 
the full (black) line is the non-spectral case, and the dashed (red) line is the maximum growth 
rated obtained when varying #0 in the spectral case. 

The essential difference between these two simulations is the number of radial modes 
that are kept in each of them. There are many more radial modes in the non-spectral case 
compared to the spectral one, in which the radial wave vector (k r ) is set by the condition of 
the field alignment of the mode Q] 



and is zero at the low field side position (poloidal angle 8 = 0). This suggests that the 
unstable modes for kgpi > 1 have a finite radial wave vector at the low field side position. It 
is well known, that a finite radial wave vector can be introduced in the ballooning transform 
through the introduction of the angle 9q such that 



The growth rate as a function of 8 for various values of kgpi is shown in Fig. |2j For 
kepi > 0.6 the most unstable mode has a finite #0 and the mode is shifted away from the low 
field side, as shown in Fig. [3j which displays the eigen function along the magnetic field for 
kepi — 1-5 and 8 = 1.2. There is no preferred sign for 80 and the mode with 8q = —1.2 is 
equally unstable but shifted in the negative 8 direction. Taking the maximum growth rate 
(by varying 8 ) for each kgpi from the spectral simulations yields the dashed (red) curve in 



k r = s8kg , 





(2) 
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Fig. [TJ There is agreement between the spectral and the non-spectral cases, as there should 




FIG. 2: (Colour on-line) Growth rates of the spectral case (7 in units Vthi/Ro) as a function of #0 
for four representative values of kgpi. The curves are denoted as follows: kgpi = 0.5 (blue) crosses, 
kgpi = 1.0 (red) circles, kgpi = 1.5 (black) pluses, kgpi = 2.0 (magenta) stars 
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FIG. 3: (Colour on-line) Eigenfunction in arbitrary units as a function of the poloidal angle 6. Full 
lines give the real, whereas dashed lines give the imaginary part. The (red) lines without symbols 
is the potential perturbation (5 (ft), the (blue) lines with the symbol V the density perturbation 
(5n), and the (magenta) lines with the closed circles is the perpendicular temperature perturbation 
(ST ± ). 



In order to get further insight into the high kg ITGs several parameter scans have been 
performed. These are done varying one parameter while keeping all the other parameters 
fixed to the standard case. The results of the R/Lt, q and s scans are shown in Figs. H] 
and |5j All calculations are performed using the non-spectral setup with kgpi = 1.9. For 
comparison, also the results for kgpi = 0.5 are shown. 

Fig. H] shows that the growth rate of the high kg ITG increases with R/Lt very similar 
to the kg pi = 0.5 mode. The mode has a higher threshold in R/Lt though is more stable 
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FIG. 4: (Colour on-line) The growth rate as a function of R/Lt (left) and the safety factor q 
(right). The (red) full line gives the result for kgpi = 1.9, while the (blue) dashed line gives the 
result for kgpi = 0.5. 

over the entire scan. In the same figure also the dependence on the safety factor is shown. 
A larger safety factor increases the connection length between the low and high field side 
and, therefore, works destabilising. This is the case for kePi = 0.5, but to a much larger 
extent for the high kg ITG. The localisation of the mode at 6 ^ makes the requirement 
of a sufficient long field line length less easy to satisfy. 




Magnetic shear Magnetic shear 



FIG. 5: (Colour ondine) The growth rate (left) and frequency (right) as a function of the magnetic 
shear. The (red) full line gives the result for kgpi = 1.9, while the (blue) dashed line gives the 
result for kgpi = 0.5. 

Fig. |5] gives the value of the growth rate as a function of the magnetic shear. The growth 
rate curve has two maxima for kgpi = 1.9, and it can be verified from the Figure of the 
frequency that these maxima belong to two different modes. A high growth rate is obtained 
only at sufficiently large shear. A high shear reduces the width of the eigenmode and is 
therefore beneficial for the # 7^ modes. The dependence of the growth rate on the inverse 
aspect ratio e (not shown) is found to be relatively weak. 



III. PHYSICAL MECHANISM 



An understanding of the physics of the high kg ITGs can be obtained by considering a 



simple fluid model. Here, the equations and normalisation given in [10j are used, and the 
reader is referred to this paper for details on the derivation. The gyro-kinetic equation, 
neglecting the parallel derivatives can be written in the form (see Eq. (68) of Ref. 10|): 

ZeV(0) 



f + v D .V/ 



v E ■ V P F M - v D ■ 



T 



M 



(3) 



where vp is the drift due to the magnetic field inhomogeneity, is the perturbed ExB 
velocity, / the perturbed distribution, the perturbed electrostatic potential, and V p is 
defined through Eq. (69) of Ref. 10|. In comparison to Ref. 10] the plasma rotation will 
be neglected, but it will not be assumed that the mode is localised on the low field side. 
Assuming a concentric circular magnetic equilibrium, 



vd ■ V = if oke cos 9 + wdK sin 6 = ikgVoK 



(4) 



where 9 is the poloidal angle, and kg (k r ) is the poloidal (radial) wave vector. In the equation 
above K is introduced to shorten the mathematics. Using Eq. (J2]) one obtains 



K = cosd + s(d-dv)sm6 
K measures the dependence of the convective derivative (vd • V) on the poloidal angle. 



(5) 



Starting from Eq. ([3]) one can follow the same procedure as outlined in Ref. [10] to obtain 
the equations for the perturbed density normalised to the background density (n), and 
perturbed temperature normalised to the background temperature (T). For singly charged 
ions, neglecting the plasma rotation the expressions are 

R „„\ 



(7) 



un + 2Kfi + 2KT = ((f)) ( 2K 



4 14 ~ / R 4 

ujT + -Kn + - KT = (0) I jL - -K 



where u is the frequency normalised to the drift frequency ujd = —kgT/eBR, and (f) here is 
the perturbed electrostatic potential normalised with e/T. Note that all terms that are due 
to the drift are proportional to K. We will therefore refer to u* D = u^K as the effective 
drift frequency. 
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The angle brackets in the equation above denote the gyro-average, or FLR effects, which 
will be modelled using a Pade approximation 

= Tuhm = Fl (8) 

where F has been introduced to shorten the notation. Finally, the gyro-kinetic Poisson 
equation is solved assuming adiabatic electrons 

n = (l + \klpi\ = G0 (9) 

where the term proportional to k\ is due to the polarization. 
From the equations above, a dispersion relation can be derived 

2 + 4 +c - (1 o) 

where 

A = G (11) 

C = ?V + G)+2 £A_M™ (13) 

3 K L T 3 K Ln 

The growth rate normalised to |a;o| (7) can be readily calculated 



V2-FT R T , R 

7 = -c^vr~ K L — (14) 

Lr V L'T ^Tcrit 



where the critical gradient is given by 



B W-^ + G^ + H^- (15) 



L<Tcrit 8 3 3 K L7V 

and we have used FG = 1. 

Fig. [6] shows the results of the growth rate, normalised to \u>d\, of the fluid model as 
a function of 9 for three values of #0 = 0, 0.5, 1.0. The left panel shows the results for 
kepi = 0.5, whereas the right panel shows the results for kgpi = 1.3. It can be seen that for 
kePi = 0.5 the mode has a maximum growth rate for 8 = 0, whereas for kepi = 1.3 the mode 
with 6*o = is stable and the most unstable mode occurs for 9 = 1.0. 

The figure of the growth rate shows that the largest growth rate at high kgpi is obtained 
for 9 ~ #o ; i-e. for k r w 0. The small value of k r then does not increase the FLR and 



FIG. 6: (Colour on-line) Growth rates of the fluid model as a function of 9 for 9q = Q (dashed), 
9q = 0.5 (solid) and 9q = 10 (dash-dotted). The left panel gives the growth rates for kgpi = 0.5, 
whereas the right panel gives the results for kgpi = 1.3. 



polarization stabilisation of the mode. Next, we clarify why the mode is strongly stabilised 
for #o = and has its maximum growth rate for 8q ^ 0. For 9 pa 6q, k r = and K = cosO. 
Therefore K — 1 at the low field side position (9 = 0) and decreases for 8^0. Eq. ( 1T^|) 
gives the dependence of the growth rate on K. If K is treated as a free parameter, and 
the density gradient is chosen to be zero for simplicity R/L^ = 0, then a maximum in the 
growth rate is obtained for 

Km = 1 R / L t (16) 

2 R/ LiTcrit 

i.e. when R/Lt > 2R/ Ltctu a maximum growth rate is obtained for the low field side 
position whereas for R/Lt < 2R/ Ltctu the maximum growth rate will be obtained for 
6 ^ 0. As kgpi is increased for fixed R/Lt, R/Ltctu increases, Km decreases, and the mode 
shifts away from the low field side. The dependence of 7 on K is shown in Fig. [7] for various 
values of kgpi. It can be seen that for kgpi = 0.5 the low field side position is the position for 
which the maximum is reached, while it is shifted away from the low field side for k$pi > 1. 

The physical reason for a maximum in K can be understood as follows. The ITG gen- 
erates ion temperature perturbations due to the perturbed ExB velocity in the background 
gradients (the term proportional to R/Lt on the right hand side of Eq. Since the drift 
(v.d) is a function of the particle energy, the temperature perturbations generate density 
perturbations through the convection (the term 2KT on the left hand side of Eq. (j6]). These 
ion density perturbations then lead to the generation of the electric field (Eq. ()9]) which is 
responsible for the perturbed ExB velocity. For K = 0, the convection due to the drift is 
zero and the mode is stable. One might therefore expect that a higher K leads to a more 
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unstable mode, and to some extent this is indeed the case, as is clear from Eq. ( 1T%|) which 
predicts 7 oc \[K. However, the Eqs. f )6|7|) also contain terms that have a stabilising effect: 
The change in kinetic energy of the ions due to the drift motion in the perturbed potential 
(the term — 4JC/3(<p) on the right hand side of Eq. (j7]), the temperature perturbations that 
are generated by the perturbed density perturbations (the term AKn/3 on the left hand 
side of Eq. (j7|)), and the fact that density and temperature perturbations have a tendency 
to propagate with different phase velocities. These stabilising terms are responsible for the 
threshold of the mode, and are all proportional to K. When the threshold is increased by 
FLR and polarization effects, and is close to R/Lt, the largest growth rate is obtained for 
K < 1, i.e. a mode shifted away from the low field side. 

,k B p = 0.9 ^ 
k e p = 2.1 

■V \ 
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FIG. 7: (Colour on-line) Left panel: Growth rates (normalised to \cjd\) of the fluid model as a 
function of K for kepi = 0.5 full (blue) line kepi = 1.3 dashed (red) line, and kgpi = 2.1 dash-dotted 
black line. Right panel: Growth rates (normalised to Vthi/Ro of spectral gyro-kinetic simulations 
with 6q = 0, with the drift multiplied with a constant a. The values of kgpi as well as the curve 
labels are the same as that of the fluid model in the left panel. 

The fluid model is, of course, a strong simplification compared with the full gyro-kinetic 
model. The fluid model not only suggests that all instabilities close to the threshold would 
have their maximum growth rate away from the low field side, it also finds no threshold for 
the ITG, since for any finite R/Lt, K can be chosen small enough that an instability arises. 
In particular the parallel dynamics (Landau damping) contained in the gyro-kinetic model 
must be considered. This stabilising mechanism is independent of K and can be expected 
to stabilise any instability for which K is too small. Nevertheless, if the explanation based 
on the fluid model is correct, its predictions should be qualitatively reproducible by the 
gyro-kinetic simulations. We discuss two tests below. 
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First, we can artificially multiply the drift velocity with a factor a (0 < a < 1), in 
spectral simulations with 6q = 0. This reduces the drift frequency and is as if we introduce 
the factor K of the fluid model into the gyro-kinetic simulations (with a = K). For those 
modes that have a maximum growth rate when the mode is shifted away from the low field 
side, one expects the maximum growth rate for 6q = to be obtained for a < 1, if the 
physics mechanism discussed above is correct. The right panel of Fig. [7] shows the growth 
rates of the gyro- kinetic simulations as a function of a for the same values of kgpi as the fluid 
model (shown in the left panel). Indeed, the gyro-kinetic simulations at high kgpi are stable 
for a = 1 and have a maximum in the growth rate for a < 1, qualitatively reproducing the 
fluid model. 

Second, as discussed above, the mechanism is not limited to kgpi > 1. Close to the 
threshold, the most unstable mode can be expected to be shifted away from the low field 
side (provided the Landau damping is small enough). Fig. [8] shows the growth rate as 
a function of 6q of the Waltz standard case with kgpi = 0.5 for several values of R/Lt 
close to the threshold of the mode. Although the effect is small, the largest growth rate is 
obtained for 8q 7^ 0. In fact, for R/Lt = 3.7, an unstable mode exists for 9q ^ whereas 
the mode at 9q = is stable, i.e. a mode shifted away from the low field side exists for 
a temperature gradient length below the threshold of the ITG obtained for 9q = 0. Both 
tests give confidence that the physical mechanism found through the analytic fluid model is 
indeed the reason for the observed behaviour of the gyro-kinetic simulations. 
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FIG. 8: (Colour on-line) Growth rates as a function of #o for the Waltz standard case with kgpi = 
0.5, and R/Lt = 3.7 open (blue) circles, R/Lt = 3.9 (black) 'x', R/Lt = 4.0 closed (magenta) 
circles, R/Lt = 4.5 (green) stars, and R/Lt = 5.0 (red) '+'. 
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IV. COMPARISON WITH PREVIOUS WORK 



Ion temperature gradient instability at sub-Larmor radius scales have previously been 
reported in the literature 2|-|6|. These modes have been found in slab geometry, as well as in 
the case of weak toroidicity. The latter condition translates to a density gradient R/L^ > 6 
for the instability to occur jjj, Such a high gradient is not usually obtained in Tokamak 
plasmas under normal operation. In contrast the high kg pi ITG described in this paper 
occurs for a wide range of R/ L n as shown in Fig. [9j and is unstable for R/L n = 0. 




FIG. 9: (Colour on-line) Growth rates as a function of RjLjq for various values of kgpi. 



There are similarities between previous work and ours. In Refs. {2, 3] it is stressed that 
the non adiabatic response of the ions at kgpi > 1 is essential for the instability to occur. 
A similar statement can be made for the modes discussed in this paper. However, the 
essential ingredient discussed in this paper, the shift of the mode away from the low field 
side, reducing the effective drift frequency, is a distinct mechanism from_that_of the works 
published to date. In particular, an inspection of the equations in Refs. 
all these references assume 8q = 0. 



3|, [5] shows that 



V. CONCLUSION 



In this paper we have shown that 

• The ITG with adiabatic electrons for standard parameters can be unstable for kgpi 
substantially larger than one. 

• Essential for this instability is a reduction of the effective drift frequency through the 
shift of the mode away from the low field side. 
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• An enhancement of the growth rate through the reduction of the effective drift fre- 
quency can be important for kePi < 1, in particular close to the threshold. 

• Unstable modes with 9 ^ can exist for ion temperature gradient lengths below the 
threshold of the mode obtained with 9 = 0. 

The existence of these modes might set additional requirements on resolution in nonlinear 

runs, and might play a role in small scale zonal flow generation. 
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